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Can a traveling wave connect two unstable states? 
The case of the nonlocal Fisher equation 
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Abstract 



This note investigates the properties of the traveling waves solutions 
of the nonlocal Fisher equation. The existence of such solutions has been 
' proved recently in [5] but their asymptotic behavior was still unclear. We 

use here a new numerical approximation of these traveling waves which 
shows that some traveling waves connect the two homogeneous steady 
^yy ' states and 1, which is a striking fact since is dynamically unstable and 

, 1 is unstable in the sense of Turing. 
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1 Introduction 

For the semilinear equation 

!^: d tU -d xxU =f(u), /(o) = /(i) = o, (i) 

a traveling wave can connect the two steady states u = and u = 1 in various sit- 
uations. This is the case when is unstable and 1 is stable (Fisher/monostable). 
This is also the case when and 1 are both stable (Allen-Cahn/bistable). It is 
known that these waves are attractive and are obtained as the long time limit 
of the dynamics ([T]) associated with compactly supported initial data (see [6]). 
When and 1 are unstable, in between there is a stable steady state of / which 
prevents any traveling wave to exists. 

For systems and non-local equations, the classification becomes more com- 
plicated because a steady state can be Turing unstable, which means that 1 
is unstable with respect to some periodic perturbations in a bounded range of 
periods. In this note we consider the nonlocal Fisher equation ([2"]). the simplest 
to produce Turing instability. 



d t u = d xx u + (m(l - (f>*u), <p*u(x)= / u{x ~ y)4>(y)dy, (2) 
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Figure 1: Numerical simulations of the time evolution for the nonlocal 
Fisher/monostable equation @ with kernel ([5]). The computational domain 
is [—80,80]. Left: the isovalues show that it is not a traveling wave but a 
more complicated structure; the bottom subplots are zooms of the top sub- 
plots. Right: the function v L connects the (dynamically) unstable state with 
a periodic tail at x — — oo. Two values of /x are used a) fi = 100; b) // = 200. 



with 

H > 0, (f>(x) > 0, (f>(0) > 0, / (f>(x) dx = 1. (3) 

Jn 

The steady state 1 can be Turing unstable when the Fourier transform of <f> 
changes sign and \i is large enough (see [SI HI [2] )- This creates several differences 
with the monostable or bistable equations. 

First of all, one can observe that the solution of ([2]) associated with a com- 
pactly supported initial datum does not converge, for large times, toward the 
traveling wave but it converges towards a more complicated structure (see [3J, 
[3] and [S]). The numerical simulation presented in b) of Figure @] shows that 
the solution of the evolution equation converges to a pulsating front, that is a 
function u(x — at,x) which is periodic in its second variable. These types of 
fronts typically arise in the framework of reaction-diffusion equations with pe- 
riodic coefficients (pQ, [7]). In a) of Figure 01 the solution seems to converge to 
the superposition of a traveling wave and a pulsating front, with two different 
speeds. These periodic patterns are a symptom of Turing instability on the full 
line E. 

Secondly, it is proved in [5] that a traveling wave u(x, t) = v(x — at) always 
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exists for all a > with a generalized formulation 

( ad x v + d xx v + fiv(l ~ cf>*v) — 0, i£l, 
} liminf ;r _ s ._ 00 v(x) > 0, v(+oo) = 0. 

The authors were only able to obtain, due to the nonlocal effect, the weak 
boundary condition at x — — oo rather than the expected condition u(— oo) = 1. 
When fi, is small enough or when the Fourier transform of the kernel F{<fi) is 
positive everywhere, the traveling wave connects to 1. But this leaves open 
the question to know whether for fj, large, the traveling wave solution of ([5]) 
connects the (dynamically) unstable state to a stable periodic state or to the 
Turing unstable state 1. Also, when F{4>) can take negative values, it is proved 
in [2] that for /x large enough monotonic traveling waves cannot exist. 
In the sequel, we will firstly focus on the case a = and 

4>{x) = \ { _ hl] {x). (5) 
Then we also test two other cases 

<K*)= (-z 2 + /3 2 )l W] , /9 = (|)* (6a) 

<P{x) = {\ + \\x\)l [ ^ 1] {x) (6b) 

In all these cases, takes negative values. For instance, for JS]), .F(</>)(£) = 

sin£/£. Finally, to complete the tests, the traveling waves for the Gaussian 
kernel 

0(a;) ^ exp (-y)' (7) 

which has positive Fourier transform, are displayed. 

In section [21 we develop a specific algorithm which allows us to build the 
traveling wave Q and not the pulsating front. Then, in section in[31 we perform 
numerical simulations to decide which alternative holds true. 



2 The algorithm 

Numerically one can only solve the problem on a bounded domain of length 
L = x r — xi (this is also the analytical construction in [5]) 

J <r L d x v L +d xx v L +fiv L {l - 4>-kv L ) = 0, xi < x < x r , 

\ v L ( Xl )=l, v L (x r )=0, v L (0) = e, {S) 

The convolution is computed by extending v L by 1 on (— oo,x;) and on 
(.T r ,+oo). The parameter e, small enough, is needed for technical reasons but 
intuitively its value has to be below the oscillations observed in Figure [5] 

Our algorithm for solving is to divide the computational domain into two 
parts: I\ — [xi,0] and 1% — [0,x r ]. Being given tr, in each interval an elliptic 
equation with Dirichlet boundary conditions is solved 

-ad x Vi = d xx Vi +fivi(l - <f>*v), Uj(0) = e, i = 1, 2, , . 
vi(xi) = 1, v 2 {x r ) =0. 
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The convolution term is computed by defining v as v\ on I\, v 2 and I 2 , 1 on 
(—00, xi) and on (x r , +00). 

But the equation does not necessarily hold true at x = 0. We define a L so 
as to impose that the jump of derivatives at zero vanishes 

a L = [d x v 2 {x r ) - d x Vi(xi)] + / fj,v(l - cj>~kv)dx. (10) 

Lemma 2.1. When {p L , , ) satisfies and U0\) simultaneously, then v L 
is C 1 on (xi,x r ) and satisfies 

We can write abstractly this problem as a fixed point for a system of two 
equations (a, v) = (Cp(u), T(cr)). It is straightforward to make it discrete using 
finite differences. The most efficient way to solve it is to use Newton iterations. 



3 The numerical results 

3.1 Convergence of the scheme 

The numerical results we present in this section are obtained with the hat func- 
tion <f> in ([5]), e = 0.1 and we study the effect of the bifurcation parameter fi. 
The diffusion term is treated implicitly by centered three point finite difference 
while the reaction term is put explicit. 

In order to verify that the iterative scheme described in section [5] converges 
to the right solution, a crucial quantity to look at is the truncation errors at 
zero 



E(0) = 



Ax 2 2Ax 

Here v^ 1; Vq , v\ are the values of v L at the grid points — Ax, 0, Ax. The conver- 
gence results are displayed in Table [1] One can see that £"(0) converges to zero 
as Ax — > 0, which shows that our numerical results is a good approximation 
of ([5]) on the whole computational domain. Note that better accuracy of E(0) 
can be obtained if we use higher order numerical integration methods for the 
convolution term. 
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20 


0.04 


3.0619 


15.3670 
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0.04 


3.0617 


15.3670 
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0.04 


3.0617 


15.3670 


20 


0.02 


1.7883 


15.4930 


40 


0.02 


1.7905 


15.5028 


80 


0.02 


1.7888 


15.4930 


20 


0.01 


0.9721 


15.5368 


40 


0.01 


0.9716 


15.5319 


80 


0.01 


0.9715 


15.5319 


20 


0.005 


0.5075 


15.5429 


40 


0.005 


0.5075 


15.5429 


80 


0.005 


0.5075 


15.5429 



Table 1: Convergence of the truncation error at zero E(0) and of the traveling 
velocity a L for various values of Ax and L, with kernel ([5]) and /i = 64. For 
L = 00 the speed is a°° — 2y/JI = 16. 
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Figure 2: The traveling wave solution for the nonlocal Fisher equation Q with 
kernel ([5]). Left: the numerical results for /i = 1 and /i = 1000. Right: the 
results for [i — 2500, the top subplot depicts v while the bottom subplot is a 
zoom of the tail. 



3.2 Convergence of the traveling waves to 1 

The traveling wave shapes for \i = 10, // = 1000 and [i — 2500 for kernel §5§ 
are depicted in Figure [5] When /i = 10, we observe a monotone traveling wave 
which connects to 1, as for the local Fisher equation. When [i grows, some 
oscillations appear. Numerically, when we increase L, the amplitudes of the tail 
decrease and the bigger fj, is, the slower the amplitudes decrease. The shapes of 
v suggest that though 1 is Turing unstable with the kernel (JSJ when [i = 1000 
and fi — 2500, the traveling waves will still connect 1 to 0. 

We do not obtain the same type of structure than when we compute the 
solution of the evolution equation depicted in Figure 01 This means that there 
exist some traveling waves that connect to 1, but that these waves do not 
attract the solution of the Cauchy problem associated with compactly supported 
initial data. 



3.3 Monotonicity of the traveling waves 

Lastly, we consider the critical value of /i for which the monotonicity of the 
traveling waves is broken. Since the monotone traveling waves always connect 
1 to 0, we perform the linearization, close to x = — oo, by assuming v ~ 1 — e Xx 
with A a real positive numberQ. After inserting this form into ((4]), using a = 2y/Ji 
and the smallness of e Xx , A can be determined by /j, through the equation 

% -2A0I-A 2 = O, (11) 



a quadratic equation for y/JI, which gives yfji = x +A ^f,^ smhA , Thus the 
critical (i that makes A no longer exist is [i c = sup A>0 ( A2+W s ^ A Asinh A ) 2 ~ 8.9. 

We have checked that this threshold /i c w 8.9 is correct on the numerical 
values. The maximum of v L is 1 when fi = 9, but exceeds slightly 1 when 

x We do not know if such a linearization is legitimate, but the technical arguments used in 
[2] to prove non-monotonicity for large /i are close to a linearization 
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Figure 3: The traveling wave solution for the nonlocal Fisher equation (jU with 
fi = 1000 and kernels in (|6]). Left: the numerical results with kernel (l6aj) . Right: 
the results for (l6b|) . In all these pictures, the top subplot depicts v while the 
bottom subplot is a zoom of the tail. 




Figure 4: The shapes of traveling wave solution for the nonlocal Fisher/KPP 
equation with Gaussian kernel Q when [i = 1000. 



/i = 10. In Figure [5] with fj, — 10, the wave is nearly monotonic, but, checking 
the numerical values, the maximum of v is 1.0028. This indicates that actually 
v might be not monotone even for L finite. 

In [5] the authors have proved that, when fi > fj, c , the traveling wave is not 
monotone. One open question is to know if the traveling wave is monotone when 
/i < \i c . Our simulation answers positively to this open question numerically. 

3.4 Traveling wave shapes for kernels in ((6j) 

The numerical results with kernels in ^ are depicted in Figure [3J With large 
/i = 1000, we can see similar phenomena such that some oscillations appear and 
the waves become nonmonotone. For the two kernels in (j6]), the state v = 1 is 
Turing unstable as for ^ and Figure [3] suggests again that the traveling waves 
will connect 1 to 0. 

Finally, to complete the tests, we show the traveling wave for the Gaussian 
kernel for fi = 1000 in Figure GO 
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